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Output regulation control for a CSTR benchmark problem is considered using a feedback linearization technique, where 
a linear control method is applied to the system for the purpose of maximizing the yield of a desired product at a specific 
operational temperature. Simulation results showed that the proposed feedback linearization-based controller strategy 
was successful in maintaining the desired product concentration at its set points, while maintaining the cooling jacket 
temperate fixed at all times, and the manipulated variables were maintained within their respective operational limits. The 
proposed feedback linearization-based controller provided very promising results, where it guaranteed a precise 
operation of the reactor with good performance in terms of a stable transition with no overshoot, and exhibited 
robustness by rejecting the tested disturbance in the form of a sinusoidal time variation in the reactant feed 
concentration. 
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Introduction 


Continuously Stirred Tank Reactors (CSTRs) are processing units that are commonly used in various continuous production, 
reaction-based, operations. Process related conditions and/or economic considerations require the operation of such reactors to be 
at an optimal point, where a desired product yield is maximized or the production cost is minimized (Edgar et al., 2001). This feat 
requires maintaining a specific set of operational conditions that dictate the processes inputs and outputs (i.e. input-output 
constraints) is undertaken by process control. One of the most widely adopted benchmark problem for the aforementioned 
situation is the CSTR reactor model presented by Klatt and Engell (Klatt and Engell, 1993; Kroll and Schulte, 2014), which 
include the reaction scheme originally presented by van de Vusse (Van de Vusse, 1964). The Klatt and Engell reactor model 
includes a realistic process that includes temperature-sensitive, consecutive and parallel liquid phase reactions, under actual 
physical considerations (Kroll and Schulte, 2014). 

In-process control, a CSTR represents an open, non-linear, reaction system model that has balance equations that combine 
kinetics and hydrodynamics. A variety of approaches are proposed for the control of the aforementioned benchmark problem in 
the literature: Chen et al (Chen et al., 1995) used nonlinear predictive control to maintain the operation of the CSTR at a desired 
set point by considering a multi-input and single-input control problems. Klatt and Engell (Klatt and Engell, 1998) have used the 
principle of gain-scheduling, while Rothfuss et al (Rothfuss et al., 1996) used a flatness-based control approach. Graichen et al 
(Graichen et al., 2004) and Perez et al (Perez et al., 2002) introduced an inversion-based feedforward control design, while 
Kvasnica et al (Kvasnica et al., 2010) proposed a method to approximate the nonlinear behavior of the system by several local 
linear models, thus allowing for a piecewise affine (PWA) model representation that predicted and optimized the reactor 
behavior. Abdalla and Shaqarin ( Abdalla and Shaqarin et al., 2017) have successfully implemented an alternative LMI (Linear 
matrix Inequality) approach on a CSTR for temperature and level control using LPV (linear parameter Varying) controller. In this 
work, state regulation control for the Klatt-Engell CSTR benchmark problem is considered using a feedback linearization 
technique, where a linear control method is applied to the system for the purpose of controlling the yield of a desired product at a 
specific operational temperature. Controlled and uncontrolled CSTR performance are compared at nominal conditions, and with 
an added disturbance in the reactant feed concentration. This alternative approach for state regulation control maintain the effect 
of the nonlinear part on the system’s performance by using exact state transformations and feedback, rather than by linear 
approximations of the dynamics. 
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1 Materials and Methods 
1.1 The CSTR benchmark model 


The benchmark model represents a constant volume CSTR with a constant volume. The reactions (Chen et al., 1995; Rothfuss et 
al., 1996) are taking place in the liquid phase, and represent the process of substance B production. This chain of irreversible 
reactions represents the production of Cyclopentenol (B) from Cyclopentadiene (A) via an acid-catalyzed electrophilic addition 
of water in dilute solution. Additionally, Dicyclopentadiene (D) is generated as a side product and Cyclopentanediol (C) as a 
consecutive product. Figure 1 illustrates the CSTR reactor in the benchmark model, along with the variables used in the 
benchmark model. 


A 
k k > A 
A>B>C (1) Oin Chin Tin re 


k2 


2A>D (2) 


where k; (1/hr) and kz (mol/m°) are reaction rate constants for 
reactions | and 2. Reactant A is used to produce B, while undesired Ou 
reactions take place producing two by-products (C and D) (Rothfuss et 
al., 1996). Both steps of reaction 1 follow the same mechanism, thus wt CAC BsT 


the same reaction rate constant (k;) is used for both steps (Klatt and T 


Engell, 1993). As shown in the figure, the continuously mixed reactor, 
with an active volume of V, is fed by reactant A at a feed flow rate of A.B, CD 
Qin (L/hr), a concentration of Chin (mol/m*), and a temperature of T; 


(°C). The feed temperature is subject to variation. None of the Fig. 1 Representation of the CSTR problem deferent variables. 
products are present in the feed. The reactor effluent is at 
temperature T (°C), with A and B concentrations of c4 and cg (mol/m*), respectively. Heat is removed by a rate of Qy (kJ/hr) to 
maintain the exothermic reactions through the usage of a cooling jacket, which has a temperature of Te (°C), and an external heat 
exchanger. 

Reactions rates for A and B are: 


ra = —kıc; — kzc2 (3) 
Tg = ky (C4 — Cp) (4) 


The dynamic behavior of the reactor can be represented by the material for A and B: 


dC,/dt = C2) (Cain = c4) — kica — kzcå 5) 
dCy/dt = -cg + k (c4 — ce) © 


In addition to the heat balance between the reactor and the cooling jacket: 


Qin 1 ARK) 
dT/dt = en Tin — T) - Cop) ka AF nas T kıcgAHRrgc T k,cxAHp 4p) + (To — T) 0) 
dTe jät = 2 y W (PT) (8) 
MxKCpy MxKCpy 


where Cp is the heat capacity of the liquid phase of the reactor, p is the density of the liquid phase, AH ,, is enthalpy of the A >B 
reaction, AHp gc is enthalpy of the B -C reaction, AH z apis enthalpy of the A >D reaction, Ag is the surface area of the cooling 
jacket, Ky is the heat transfer coefficient for the cooling jacket, mg is the coolant mass, and Cp;is the heat capacity of the coolant. 
The reaction rates in Equations 3 and 4 are affected by temperature through the Arrhenius law: 
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k; = k? exp ( = ) (9) 


T+273.5 


where i denotes reaction rate constants 1 and 2, k? is the collision factor for rate constant i, and Æ; is the activation energy for 
each of the aforementioned reaction constants (Chen et al., 1995). Concentrations of substances C and D (the by-products) are 
not considered in the context of the model, and they do not influence the dynamics of the process (Rothfuss et al., 1996). Table 1 
illustrates the values for the different physical and chemical parameters in the benchmark model. 


As it is shown in Table 1, values of the 
collision factors and Enthalpy’s do indeed 


Table 1 Values for the different physical and chemical parameters in the benchmark model (Chen 


et al., 1995) 


have some uncertainty. This work used Symbol Parameter Value 
values of the best estimate in the v Reactor Volume 0.01 m? 
simulations. c, Heat capacity of the liquid phase of the 3.01 KJ/kg.K 
reactor 
1.2 The control problem p Density of the liquid phase 934.2 kg/m? 
AHR 4B Enthalpy of the A -B reaction 4.2 +2.36 kJ/mol A 
The control problem of the specific AHR gc Enthalpy of the B >C reaction -11.0 + 1.92 kJ/mol B 
benchmark model CSTR in hand involves AHrap Enthalpy of the A >D reaction -41.85 + 1.41 kJ/mol B 
a constrained control problem of a Multi- Ar Surface area of the cooling jacket 0.215 m? 
input Multi-output (MIMO) coupled Ky nee transfer coefficient for the cooling — 4032 kJ/hr. K.m? 
nonlinear system. Feedback linearization Jacket 
technique is a good candidate for such ”* Coolinpmass Toke 
problem; since it is capable of canceling C» Heat capacity of the coolant 2.0 kJ/kg.K 
the nonlinearities (fully or partially) and/or k? Collesion factor for rate conatant 1 1.287 + 0.04 x 10”? 1/hr 
6 
to decouple the model equations. A state k? Collesion factor for rate conatant 2 NRS vas x 10 
feedback controller var be designed based Ei Activation energy for rate conatant 1 -9758.3 K 
on the feedback linearized system. saan 
E Activation energy for rate conatant 2 -8560.0 K 


1.3 Feedback linearization 


Feedback linearization is a technique used for nonlinear systems control design. Basically, the nonlinear dynamics of the system 
to be controlled is algebraically transformed, fully, or partially into a linear one. This enables the usage of linear control 
techniques in the system’s controller (Slotine and Li, 1991). The value of this approach, compared to conventional linearization, 
stems from the fact that feedback linearization is achieved by exact state transformations and feedback, rather than by linear 
approximations of the dynamics. Thus, maintaining the effect of the nonlinear part on the system’s performance. The mechanistic 
approach used in feedback linearization, i.e., canceling the nonlinearities and imposing a desired linear dynamics, can be simply 
applied to a class of nonlinear systems described by the so-called companion form. A system is said to be in companion form if 
its dynamics is represented by: 
x™ = f(x) + b(x)u (10) 

where u is the scalar control input, x is states, x = [x, x, ..., x’ is the state vector, and f(x) and b(x) are nonlinear functions of 
the states. This form of a dynamic system model is unique in the fact that it is easily transformed into a linearized system via 
direct feedback linearization. 


The companion from Equation 10 can be transformed to a state-space representation as follows. 


xX, X2 

X| _ X3 (11) 
x, (x) + b(x)u 
The nonlinear state-space is easily transformed to linear one by selecting an appropriate control law in the form of: 
u=-[v-f] (12) 


The closed-loop system now is transformed to an equivalent linear state-space equation as follows: 
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x" =v (13) 


where v is the linear control law that can be designed via, state-feedback control, optimal control, pole-placement, ... etc. The 
result is a linear input-output map between the new input and the output. Figure 2 illustrates the input state linearization general 
form. 

The resulting system structure is unique and system-specific, and most of the process systems cannot be represented as such, for 
instance, the CSTR. Albeit, feedback linearization can be used to fully (Shagarin and Abdalla, 2009) or partially (Shaqarin et al., 
2014) cancel the nonlinearity for CSTR model, and remove the coupling between the state variables, especially, the controlled 
variables. 

In the CSTR problem case, the system state is x= [c+ cp, T, T,]. Controlled 
state variables are the effluent substance B concentration (cg) and the jacket ' ) uul 
temperature (Tc) (Rothfuss et al., 1996). The cooling jacket temperature is 
controlled as a means to control the reactor temperature. The manipulated 
variables are the reactor volume normalized flow rate Q;,/V (1/hr), and the 
heat removal rate by the jacket of Q, (kJ/hr). From an operational point, 


Linearization loop 


Qin/V is manipulated by changing the flow rate into and from the reactor, 
while Qy can be manipulated through changes in the coolant flow rate and/or 
using a multiple heat exchangers setup, to name some. However, these changes 


can be strictly performed within specific constraints dictated by the process 
itself. The constraints for the Klatt-Engell CSTR benchmark model 
manipulated parameters are shown in the next section. 


Fig. 2 General form for the input state linearization 


1.4 CSTR operation 


The state regulation problem requires operating the reactor at set points for Table 2 Constraints on the state variables and 
manipulated variables (Kvasnica et al., 


substance B effluent concentration (cg), while maintaining the temperature of the 2010; Kroll and Schulte, 2014) 
cooling jacket (Tc) at a desired value. Table 2 shows the constraints on the state 
variables (for safe operation) and manipulated variables (Kvasnica et al., 2010; Kroll Model variable Constraints 
and Schulte, 2014). In this work, the setpoints for cz is 1.09 x 10° (mol/m*), which is Qin/V 5 to 35 
considered as the point where the maximum yield of B is obtained (Chen et al, o, -8500 to 0 
1995). Values for the remaining model variables required to achieve this point (max g 2.6 x 10° to 3.7 x 
yield) are shown in Table 3. The cooling jacket temperature is also regulated at 109 á 10° (mol/m’°) 
(°C), which is lower than the temperature shown in Table 2. Finally, the CSTR CR Hen a to 1.1 x 
model performance was simulated using MATLAB (The MathWorks Inc., Natick, T momy 
MA). 80 to 150 °C) 
Tc 

2 Results and Discussion Table 3 Operational parameters at maximum, 
21 Open-loop operation: optimal, yield for B (Chen et al., 1995). 

Model inputs Optimal value 
In order to validate the performance of the CSTR model, a simulation was Cain 5.1 x 10° (mol/m*) 
performed for the maximum B yield point mentioned in Table 2. Figure 3 illustrates Tn 104.9 °C 
four plots that represent the simulation result for an open, uncontrolled CSTR. @Q,,/V 14.19 (1/hr) 
Initial values for the model state variables in the reactor were assumed to be zero g, -1113.5 (kJ/hr) 
(Ca, Cg, T, and Tc= 0). After 2 hrs of simulated operation, ueni steady-state G 2.14 x 10° (mol/m?) 
concentrations of A and B were 2.113x 10° and 1.09 x 10° mol/m’, respectively, i 1.09 x 10° (mol/m’) 
corresponding to a yield of 0.516 for B. Reactor and jacket steady-state T 142°C 
temperatures were 114.25 and 112.93 °C, respectively. These results were identical Te 112916 


to the values reported by Chen et al., (1995) and Rothfuss et al., (1996). 
In order to gain a better understating regarding the performance of the benchmark CSTR problem, multiple simulations were 
performed using values of the model’s manipulated variables within the range depicted in Table 1. Figure 4 shows four plots 
representing the resulted contour plot for the outcome of the simulation. The figure depicts how each of the state variables varies 
within the permissible range of the manipulated variables. 
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o 1 1 1 1 1 1 1 
0.0 0.3 DS 0.8 1.0 1.3 1.5 1.8 2.0 


Time (hrs) 
Fig. 3 Simulation results for the CSTR benchmark model at Cain=5.1x10° (mol/m*), Tin =104.9°C, Qn/V=14.19 (1/hr), and Qy=-1113.5 (kJ/hr), where A: 
reactant A effluent concentration, B: Substance B effluent concentration, C: effluent stream temperature, and D: cooling jacket temperature. 


Ca 
=a 
Ea 
£ li 
Š =z 
= = 
= mz 
|| 
a 
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m a: 
= = = 
s = az 
= = = 
& = = 
HE 110 Es 
Ez 


10 20 30 40 
O,,/V (/hr) O,,/V (hry 
Fig. 4 Contour plot representing the simulated operation of the CSTR benchmark model at Cain=5.1X10° (mol/m*) and T;,=104.9°C with the manipulated 
variables 5<Q;,/VS 35 (1/hr), and -8500<Qy<0 (kJ/hr). 


As seen in the figure, the maximum concentration of B results from operating with a Q,,/V between about 9 and 25hr! and a Quy 
between about -4200 and O (kJ/hr), while maintaining the constraints on the manipulated variables shown in Table 1. This range 
corresponded to a reactor temperature and a cooling jacket temperature between about 90 and 117°C, which for the cooling jacket 
temperatures encompasses the desired set point for it at a temperature of 109°C. 


2.2 Output regulation Control 


Typical simulations of the closed-loop nonlinear CSTR system response with feedback-linearization and state-feedback 
controller are discussed herein. The set-up understudy is to operate the reactor for the purpose of maximizing the yield of the 
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desired product B at a specific cooling temperature. Plots in Figure 5 (solid lines) shows the controlled response of the CSTR 
model with feedback linearization for the operation at a set-point for cg at 0.6x10° and 1.09x10° (mol/m’), and a set-point for Tc 
at 109°C, at inputs of cam=5.1x10° (mol/m*) and Tj,=104.9°C. 


4500 L 1299 
4000 


3500 1000 


W 
S 
i=) 
© 


(mol/m*) 


800 


Cp (mol/m*) 


= 2500 
2000 
1500 


c 


600 


114 
113 
112 
111 
110 
109 
108 
107 


T CO 
T, CC) 


Controlled 
Uncontrolled 


Q,/V (\/hr) 
w e 
w oO 
Qy (kJ/hr) 


0.0 0.3 0.5 0.8 1.0 1.3 1.5 1.8 2.0 0.0 0.3 0.5 0.8 1.0 1.3 1.5 1.8 2.0 
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Fig. 5 Solid lines show the simulated controlled operation of the CSTR model with feedback linearization, with the set set-point for cg at 0.6x10° and 
1.09x10° (mol/m’°), and Tc maintained at 109°C, while dashed lines shows the simulated uncontrolled operation of the CSTR at values of Q;,/V 
and Qy manually changed between 63.4 and 14.6 (1/hr) and -1090 and -2090kJ/hr, respectively; at Cam=5.1x10° (mol/m°) and T;,=104.9°C. 


For this nominal scenario, as seen in the figure, the controller was able to regulate the operation of the CSTR, with a settling time 
of ~ 0.2 (hr), to produce cz at the desired maximum yield point. This was achieved by maintaining the desired T, at the given set- 
point of 109 °C). The control effort remained within the constraints sanctioned on the system, changing Q,,/V between 63.4 and 
14.6 (1/hr) and Qy between -1090 and -2090kJ/hr. The same can be said for the remaining state variables, as they remained 
within the safe operation range (Table 1). The regulation of cg at its maximum yield point was possible with a Tc value that is 
lower than the required temperature for optimal operation (Table 2). This was possible by increasing the heat removal rate at the 
jacket (Qy), which increased from -1113.2kJ/hr, to -2090.8kJ/hr. This change slightly reduced the reactor temperature to 111.3°C. 
On the other hand, values of Q;,/V settled back to a value that is almost identical to the optimal point (Table 2). The effectiveness 
of the proposed control technique can also be shown in Figure 5 (dashed lines), where the simulated uncontrolled operation of 
the CSTR is depicted. This simulation was performed with the same input of Cain=5-1x10°mol/m’, and T;,=104.9°C. Values of 
Qin/V and Qy were manually changed to mimic the range resulting from the controlled operation to be between 63.4 and 14.6 
(1/hr) and -1090 and -2090kJ/hr, respectively. Compared to the controlled operation, the response of cg and T, reached and 
maintained their desired set points values. However, the controller allowed for a better performance in terms of overshoot and 
settling time for both outputs. 

To test the robustness of the designed controller in disturbance rejection, the regulation control problem was simulated with 
sinusoidal disturbance as time variation in the inlet concentration substance A concentration (c,;,) in the feed by +100mol/m’. 
Perturbations in c, ;, where achieved by changing its value to be: c4 ;,=5100+100 xsin (8x?) (mol/m*). Plots in Figure 6 shows the 
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simulated controlled response (solid lines) and uncontrolled response (dashed lines) for the CSTR benchmark model with 
feedback linearization in the presence of the external disturbance (c,;,). 


£ Controlled 

© 

E 
5 Uncontrolled 
g 

o 


1200 

E “E 1000 
E È 

E = 800 

600 

113 113 

112 112 

© lll © ill 

= 110 = er 110 


ji 
S 
O 
om 
S 
O 


-500 
-1000 
-1500 
-2000 
-2500 


-3000 
0.0 0.3 0.5 0.8 1.0 1.3 1.5 1.8 2.0 0.0 0.3 0.5 0.8 1.0 1.3 1.5 1.8 2.0 
Time (hrs) Time (hrs) 
Fig. 6 Solid lines show the simulated controlled operation of the CSTR model with feedback linearization, with the set set-point for cg at 0.6x10° and 


1.09x10°(mol/m? ), and Tc maintained at 109°C, while dashed lines shows the simulated uncontrolled operation of the CSTR at values of Q;,/V and 
Qu manually changed between 63.4 and 14.6 (1/hr) and -1090 and -2090kJ/hr, respectively; at c4in=5100+100sin (8xt) (mol/m*) and T;,=104.9°C. 
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in 
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Despite the external disturbance in the inlet reactant concentration in the CSTR feed, the controller was successful in regulating 
the outputs to the set-points. This observation leads to the conclusion that the controller is robust in disturbance rejection if it is 
found in the form of time variation of the inlet concentration. Compared to the controlled CSTR behavior, subjecting the 
uncontrolled CSTR to the same disturbance (dashed lines), resulted in a great degree of disturbance in the behavior of all the 
model outputs. This clearly shown in Figure 6, which shows how the outputs did not stabilize but rather contaminated with the 
same frequency of the external disturbance in cám. Predictably, the current controller improved the performance significantly and 
was more robust against external disturbance compared to open-loop control. 
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Conclusions 


This study aimed to introduce an alternative approach for controlling an of a complex, highly nonlinear, benchmark chemical 
reactor, that included a CSTR operating with temperature-sensitive, sequential and parallel reactions. Simulation results show 
that the proposed feedback linearization-based controller strategy was successful in regulating the desired product concentration 
to its set points while maintaining the cooling jacket temperate fixed at all times. The proposed feedback linearization-based 
controller provided very promising results. Not only does the controller guarantee a precise operation of the reactor, it also 
provided, compared to uncontrolled operation, improved performance in terms of settling time and overshoot. Disturbances in the 
feed reactant concentration did not affect the performance of the controller. The manipulated variables were maintained within 
the operational limit for each one of them, and the controller was robust against the tested external disturbance. 


Nomenclature 

Acronyms 

CSTR Continuous stirred-tank reactor 

Symbols 

A: =Reacting substance [-] 

B: =Intermediate, desired product [-] 

C: =Final product 1 [-] 

D: =Final product 2 [-] 

kir: =Reaction 1 rate constant [1/hr] 

kz: =Reaction 2 rate constant [m°/hr.mol] 
V: =Reactor volume [m°] 

Qin: =Feed volumetric flow rate into the reactor ({m*/hr] 

Ou: =Jacket heat removal rate [kJ/hr] 

Caine =Substance A concentration in the reactor feed [mol/m>] 
CAL =Substance A concentration in the reactor and effluent [mol/m>] 
CB: =Substance B concentration in the reactor and effluent [mol/m3]1] 
Tint =Reactor feed temperature [°C] 

T: =Reactor and effluent temperature [°C] 

Te: =Cooling jacket temperature [°C] 

Ta: =Overall reaction rate for A [mol/hr.m*] 
rp: =Overall reaction rate for B [mol/hr.m>] 
AHR ap: =Enthalpy for the A >B reaction [kJ/mol A] 
Ar ac: =Enthalpy for the A >C reaction [kJ/mol A] 
AHr gc: =Enthalpy for the B>C reaction [kJ/mol B] 
Cp: =Heat capacity of the liquid phase of the reactor [kJ/kg.K] 
Cp;: =Heat capacity of the liquid coolant [kJ/kg.K] 
Ar: =Surface area of the cooling jacket [m?] 

Kw: =The heat transfer coefficient of the cooling jacket [kJ/hr.K.m’] 
mg: =Coolant mass [kg] 

k}: =collision factor for rate constant 1 [hr] 

k}: =collision factor for rate constant 2 [m*/hr.mol] 
E;: =Reaction / activation energy [K] 

x? =State variables vector [-] 

fandb:  =Functions [-] 


Greek letters 
p: =Density of the reactor liquid phase [kg/m] 
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